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1  It  is  known  that  a  substantial  increase  in  the  load  capacity  of 
a  self-generating  gas  bearing  can  be  obtained  at  a  correspond¬ 
ing  increase  in  the  gas  pressure  on  the  bearing  ends  only. 

This  can  be  realized  structurally  by  blowing  the  gas  into  the 
bearing  trough  through  two  annular  grooves  in  the  bushing, 
by  locating  the  grooves  near  by  the  ends,  and  by  maintaining 
a  constant  and  equal  gas  pressure  in  them.  In  this  case, 
the  middle  section  of  the  bearing,  inclosed  between  the 
grooves,  will  act  as  a  self-generating  bearing,  and  the  end 
sections  will  operate  under  conditions  of  an  axial  pressure 
drop.  „A  method  of  calculating  the  Reynolds  equation  is  pre¬ 
sented,  permitting  a  solution  of  the  problem  to  be  obtained 
in  relatively  simple  close  approximations.  An  analytical 
solution  of  the  equations  of  successive  approximations  can  be 
plotted  for  the  middle  section  of  the  bearing,  but  the 
pressure  determination  in  the  end  sections  generally  requires 
the  numerical  integration  of  the  ordinary  differential  equations 
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THEORY  OF  A  RADIAL  GAS  BEARING  WITH 
CIRCULAR  SUPPLY  GROOVES 

[Article  by  A.  I.  Snopov;  Moscow,  Doklad  na  Soveshchanll  po 
Gazovov  Smazke  Podshipnikov,  Russian,  12-14  February  19&o, 

99  bT*m — 


It  is  known  that  an  essential  increase  in  the  load- 
bearing  capacity  of  a  self-generating  gas  bearing  may  be 
achieved  only  with  a  corresponding  increase  of  gas  pressure 
on  the  bearing  face.  Structurally,  this  may  be  accomplished 
if  the  gas  inblow  into  the  bearing  is  done  through  two  circu¬ 
lar  grooves  on  the  bushing,  situated  close  to  the  end  faces, 
and  if  the  pressure  of  feeding  gas  is  maintained  stable  and 
constant.  In  this  case,  the  middle  portion  of  the  bearing 
between  the  grooves  will  work  as  a  self-generating  bearing, 
and  the  end  segments  will  work  under  conditions  of  axial 
pressure  drop. 

2  2 

Given  is  a  method  of  solving  Reynolds  equation  -  p  h 
which  determines  pressure  distribution  in  the  lubricating 
layer.  This  method  permits  obtaining  a  relatively  simple 
solution  of  the  problem  in  high  approximations.  For  the 
middle  part  of  the  bearing,  an  analytical  solution  of  the 
equations  of  consecutive  approximations  may  be  structured, 
and  the  determination  of  pressures  at  the  end  sections  re¬ 
quires,  in  a  general  case,  numerical  integration  of  the 
ordinary  differential  equations  appearing  at  each  Iteration 
stage . 

1 .  Formulation  of  the  Problem  and  Basic  Equations 

Design  computation  of  a  radial  bearing  with  several 
circular  supply  grooves  requires  the  determination  of  pressure 
distribution  both  in  the  sections  between  grooves  and  at  the 
end  sections.  We  -shall  examine  one  such  section.  Let  L  be 
the  section  length,  r  the  shaft  radius,  rj  the  bearing  ra¬ 
dius,  OJ  the  angular  velocity  of  shaft  rotation,  and  p*  and 
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With  this,  equations  and  boundary  conditions  assume  the  form 

**  o.?. 


£  U9 

d$z  =  % 


d$ 


*/• 


&$*  =  *dz 
#{/>«?)  d{Paz) 
d(f  _  dZ 


*$ 

o> 


o  > 


“<>  =  / 


u/>  -  aZ  -  u9  9  ° 

/>*  /  at-  2=  0,p*X 

We  shall  Introduce  new  variables 


S'** 

a  4~S«*9, 
at  ar  -  e. 


Z*{§  s'-'S(?)y  >  ?s  ?*>  >  J>  =  /+pco$9. 


(5) 


(6) 


With  these  variables,  equations  (4)  and  boundary  conditions 
(5)  are  written  thus 


From  the  first  two  equations  (7)  taking  into  account  (8)  we 
find  that 


“ft*  rj .  <i  -  / . 

where  - 

0/7  +  &?ps;n<p ,  £=#/>&  •  (9 

We  shall  multiply  th^,  third  equation  (7)  by  and 

we  shall  integrate  it  for  ?  within  the  limits  from  0  to  1 . 
After  some  simple  computations  taking  into  account  (8)  we 
obtain  equation 
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(10) 


d<f  d<p 

which,  upon  substitution  into  it  of  the  values  of  magnitudes 
a  and  b,  we  shall  convert  to  the  form 

H*  «  oA-n  (11) 

f2?pe*r  * 

where  6  6/ucO  / I  ^ 

A  =  T  =  “a"/  f' 


2.  p2h2  -  Method 

The  form  of  equation  (11)  naturally  suggests  accepting 
q2  =  p2h-  as  the  unknown  magnitude  and  substituting  it  anci 
the  magnitude  q  in  the  form  of  series  according  to  th-  d„g  — 
of  the  relative  eccentricity 

t’iy*  >’£*''•  <12) 

From  the  condition  /  £7  •/)  c-»  -  fyn 


we  find  that 


.  -  . 

$,•  ~-js° - £,mw ■  (i 

In  agreement  with  (11),  functions  sn  are  determined 
sequentially  from  equations 


£V £'  „  )- 

diz  S(pz  d<plf0/~  fylZff*?*'*/ 

/d\.f  \  dsn.j  , 

where  functions  sn  should  satisfy  conditions 


sos*>  s/ - 2co* <? >  cos2<p,  S^Ofes)9*'  2*0, k 

Wl  s^ofos)61  z,e.  (15) 


Assuming  that  s0  =  sQ(z),  we  find  easily 

s* ' A  e§  '  frr('-  ej)‘  > 

At  n  =  1,  equation  (14)  assumes  the  form 


£%  ^  *%  a  ds, 

dZZ  d<pz  <j,0  dp 


-  2s0cos(p 


We  assume 


st  x  crf(z)sln  (p  /■  i(  (z)  cos  <p  . 


(16) 


(17) 


(18) 


On  the  basis  of  (15)  and  (17),  we  have 

af  ~  Cff  *  “*  if  *  O  £*-  g  -  —  a  -  -  Zsn 

%  f  '  °  (19) 

ai(°)  *  ^  £,{0j-2,  i/t)-2t*. 

We  shall  introduce  a  complex  function 

xSf  =  ctf  +  ii, . 

This  function,  in  agreement  with  (19)  satisfies  conditions 


(20) 


J-J 


’Jo*  Mf(o)'  o. 


(21) 


Assuming  that  in  (14)  n  =  2,  we  obtain 

dz*  dp*  °o  o<?  L  «p 

-  j  <7/-  af]sh2y  +[“2y3afi  ~2  %  ~  &f]cos2<?  * 


(22) 


If  it  is  assumed 

•  Sg  »  ctz(z)ci/>  2 <p  +  &2(z)co$2(p  +  c2  , 


then 


To  «f0 
2JI  ..  A 


(23) 


(24) 


Jo 

c2  *  -  j  ££  > 

az(0)  -  a/#«  0,  4^*/ ’ fyQ*il*>Ci(0H*c2(tJmilz. 


We  shall  introduce  the  function 

+  &-{&  ’ 


-  a2 


(25) 


evidently,  it  is  determined  in  relation  to  (24)  from  condi¬ 
tions  : 


2-f0 "*j**',*-s 


0» 

(26) 


We  find  easily  also  that 


(27) 


We  shall  limit  ourselves  by  the  indicated  two  approxi¬ 
mate  values.  Knowing  wj  and  v.'2  we  may  determine  C4  and  q2* 

As  the  result  we  find  that 

/  ^5Srf>^(28) 

1  -  - 


where 


^•on 


}<  *  JZ0S/  "  ^ f  *  » 


+  i 8.  -  —  zif.  +  q  1  ,  ft. +>  k 
'  '  2o  7  Jo  <  ' 

fo 


/ j 


>0 


<5  «  >  *?  '■zfo* 


e  : 


In  the  general  case  when  *  ,  it  will  not  be  pos¬ 
sible  to  obtain  an  analytical  solution  of  equations  (21)  and 
(26)  in  a  closed  form,  and  for  practical  application  one  may 
take  advantage  of  their  numerical  integration.  In  the  case 
/  ,  the  solution  may  be  presented  in  the  hyperbolic  func¬ 
tions  of  a  complex  argument  [2],  at  which  point  we  shall  stop. 
We  should  note  that  in  its  first  approximation  it  is  in  agree¬ 
ment  with  Osman's  solution  [l]  . 


3.  Lubrication  Effect  on  the  Shaft  and  Gas  Consumption 

The  main  vector  of  forces  of  pressure  applied  to  the 
examined  section  of  shaft  has  the  following  components 

(  e  2*  . 

(29) 

00  00 

and  we  shall  represent  it  in  the  following  complex  form 

-  (l-Uu.B, :>  4  nr.&J. 

P  J‘  H1  ’  '  ’J 

-  4  ■*  ^-^4:  **?**-£  • 


(30) 


The  moment  of  the  forces  of  friction  applied  to  the 
shaft  along  section  L,  with  respect  to  the  z  axis  may  be,  as 
usual,  represented  by  a  formula 


M 


27ir03 JL/u  co 

~  2  • 


(31) 
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Computing  the  consumption  of  lubricant  per  second  in 
the  examined  section  from  the  formula 

r,  **  ecosy 

Q  -  J  J  J0x/,gt'c//'c/(p  , 

o  ft, 

we  find  that  it  is  independent  of  the  angular  velocity  of 
shaft  rotation  and  equals 


q  A  fj,  *i*  ) 

3/ui  fi,  [**»?*-)’ 


(33) 


Since  at  A  =  0,  w^  =  W2  =  0,  then,  in  agreement  with 
(30)  and  Px  =  py  =  °»  consequently,  in  a  bearing  with  a  circu¬ 
lar  inblow,  load-bearing  capacity  is  produced  by  the  shaft 
rotation. 

4.  Case  of  High  Velocities  of  Shaft  Rotation 

In  the  case  of  high  velocities  of  shaft  rotation,  A. 
is  large  and  it  is  possible  to  construct  an  asymptotic  solution 
of  equations  (21)  and  (26).  It  will  contain  the  functions 
which  are  solutions  of  equations  at  %  s  oo  and  boundary  layer 
functions  of  the  form  g-tfC/eg  *  When  computing 

the  integral  (30),  in  the  first  approximation  one  may  neglect 
the  squares  of  boundary  layer  functions  and  retain  only  the 
integrals  of  threshold  functions  corresponding  to  ^ 

Assuming  in  (21)  and  (26)  A  a  oo,  we  easily  find 
that  the  threshold  solutions  have  the  form 


*£  2irz 


Jo  * 


“I?/ 


(34) 


At  the  same  time 


(35) 


But 


?M- iH'-»  7- 
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consequently,  in  the  case  of  high  velocity  of  the  shaft  rota¬ 
tion 


where 


while 


*(*>' -h £)*]  >  *-'-(£■/ 

X(o)  =  /  ,  Kft)  =  j  ■ 
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